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The dynamics of cold trapped ions in a high-finesse resonator results from the interplay between 
the long-range Coulomb repulsion and the cavity-induced interactions. The latter are due to multiple 
scatterings of laser photons inside the cavity and become relevant when the laser pump is sufficiently 
strong to overcome photon decay. We study the stationary states of ions coupled with a mode of a 
standing-wave cavity as a function of the cavity and laser parameters, when the typical length scales 
of the two self-organizing processes. Coulomb crystallization and photon-mediated interactions, 
are incommensurate. The dynamics are frustrated and in specific limiting cases can be cast in 
terms of the Frenkel-Kontorova model, which reproduces features of friction in one dimension. We 
numerically recover the sliding and pinned phases. For strong cavity nonlinearities, they are in 
general separated by bistable regions where superlubric and stick-slip dynamics coexist. The cavity, 
moreover, acts as a thermal reservoir and can cool the chain vibrations to temperatures controlled by 
the cavity parameters and by the ions phase. These features are imprinted in the radiation emitted 
by the cavity, which is readily measurable in state-of-art setups of cavity quantum electrodynamics. 

PACS numbers: 42.50.Ct, 37.30.+i, 63.70.+h, 61.44.Fw 


Cavity quantum electrodynamics (CQED) with cold 
atomic ensembles provides exciting settings in which to 
study the physics of long-range interacting systems mm- 
The latter are found when the interparticle potential 
in D dimensions exhibits a scaling with the distance r 
slower than 1/r^ [3]. This property is of relevance from 
the nuclear scale to astrophysical plasmas and leads to 
non-additivity of the energy, whose consequences are, 
amongst others, ensemble inequivalence and metastable 
states with diverging lifetimes [3]. 

The dynamics of atoms in single-mode high-finesse res¬ 
onators exhibits several analogies with well-known theo¬ 
retical models of the statistical mechanics of long-range 
interacting systems [4]. High-finesse cavities, in fact, 
trap photons for a sufficiently long time so that mul¬ 
tiple scatterings can occur among the atoms inside the 
resonator. This gives rise to an effective interatomic po¬ 
tential whose range can scale with the size of the system 
[5] and that can induce spontaneous atom ordering [6l |7] 
even in one dimension [8]. Additionally, the system is 
intrinsically out-of-equilibrium since the resonator dissi¬ 
pates light, thus non-trivial phases are observed only in 
the presence of an external drive. Under these premises, 
photon shot noise gives rise to global retardation effects, 
which can effectively cool the atomic motion ninH]. 

The interparticle forces in general compete with these 
dynamics. When the interactions are short-ranged, at 
ultra-low temperatures their interplay with the cavity po¬ 
tential can give rise to exotic phases, which tend to max¬ 
imize photon scattering into the resonator and thus the 
strength of the long-range intracavity potential [mail- 
ill]. It is to a large extent unknown, however, how these 
dynamics are modified as the range of the competing in¬ 


teracting potential is increased. This question acquires 
further relevance in view of experimental setups trapping 
cold ions within high-finesse cavities [HHH]. 

In this Letter we theoretically characterize the effect 
of cavity back-action in the presence of the compet¬ 
ing Coulomb interaction between N ions with the same 
charge q and mass m. The ions are confined by an exter¬ 
nal trap inside a standing-wave resonator of wavelength 
A, in the geometry illustrated in Eig. [^a), where their 
motion is assumed to be one-dimensional along the x- 
axis. This setup is expected to simulate the Erenkel- 
Kontorova (EK) model [T9H2T] which describes a chain 
of elastically-bound particles subjected to an external pe¬ 
riodic potential (substrate) in one dimension [121 HSl HU" 
|25] and reproduces the salient features of stick-slip mo¬ 
tion between two surfaces. When the periodicity A/2 
of the cavity optical lattice is incommensurate with the 
characteristic interparticle distance d of the ions (and the 
cavity nonlinearity is negligible), the ions ground state 
can be either sliding or pinned: In the sliding phase the 
forces giving rise to sticking can cancel, so that the min¬ 
imal force for initiating sliding vanishes. Static friction 
becomes significant when ions are pinned by the cavity 
potential, in this case their distribution can still be in¬ 
commensurate with the lattice periodicity, exhibiting de¬ 
fects (kinks). In the EK model the sliding and pinned 
phases are separated by the Aubry transition, whose con¬ 
trol field is the relative amplitude of the periodic poten¬ 
tial [26ll28] . and whose hallmark is the abrupt growth 
of the lowest phonon frequency (phonon gap) [29]. In a 
finite chain with free ends and an inversion-symmetric 
potential, this transition is characterized by symmetry 
breaking [2TJ|28]. 
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FIG. 1: (a) An array of cold ions in the optical lattice of 
a high-finesse cavity is an exotic realization of the Frenkel- 
Kontorova model. The cavity is pumped by a laser with am¬ 
plitude rj and decays at rate k, the intracavity photon number 
is determined by the ion density: this gives rise to a globally 
deformable potential which depends nonlinearly on the ions 
positions, (b) The functional form of the cavity-induced po¬ 
tential for one particle: In the pinned phase, for \C\ < 1 
the shape is approximately sinusoidal, for \C\ > 1 it becomes 
flat everywhere except in the vicinity of minima (C > 0 ) or 
maxima (C < 0). (c) Sketch of the phase diagram for the sta¬ 
tionary state as a function of rj and of the cavity nonlinearity 
C (cooperativity). Typically for \C\ > 1 bistable phases can 
be observed (hatched region) signifying that superlubric or 
stick-slip dynamics are found depending on the variation of 77 
with time. 


This behaviour is substantially modified when the pe¬ 
riodic potential of the resonator mediates a long-range 
interaction among the atoms in the dispersive regime of 
CQED [UlSOl. Here, the coupling at strength to a cavity 
field mode, with spatial mode function cos{kx) and k = 
27VI induces the conservative potential hUoncos^{kx) 
for sufficiently large detuning Aq = cjc — (^e\ between the 
frequencies of the cavity mode and ion dipolar transi¬ 
tion. Then, the potential amplitude is proportional to 
the intracavity photon number operator n. The strong¬ 
coupling regime of CQED is reached when the dynam¬ 
ical Stark shift per atom Uq = g‘^/^o is such that 
N\Uo\/k > 1 , with K the cavity loss rate [Tl[30]. In this 
regime, the mean intracavity-photon number n = (n), 
and thus the depth of the cavity optical lattice, is a non¬ 
linear function of the ions positions m- 

n=\7]f/{K^ + Aes{{Xj}f), ( 1 ) 

where 77 is the amplitude of the driving field and 

Aeff= Ae - NUoBnUxj}) . (2) 


The detuning Ac = ujp — uJc of the pump from the cav¬ 
ity frequency is thus shifted by NUoB]sf{{xj}), where the 
function B^ = Xlj cos‘^{kxj)/N depends on the ions po¬ 
sitions Xj along the cavity axis and is the so-called bunch¬ 
ing parameter, as it measures their localization at the 
potential minima. This frequency shift is at the origin of 
the nonlinear dependence of n on {xj} and gives rise to a 
deformation of the effective potential the ions experience, 
as shown in Eig. Bb), which can be expressed in terms 
of the effective mean-field potential 

Kav = arctan(Aeff({a;j})/K), (3) 

and whose derivation is reported in Refs. [211 Eg. Its 
functional dependence is reminiscent of the nonlinearly 
deformable potential discussed in Refs. [221 El]- Poten¬ 
tial Vcav, however, also mediates a multi-body long-range 
interaction between the ions and thus acts as a globally 
deformable potential, since the potential depth depends 
on the global variable B^- It competes with the poten¬ 
tial Eion due to the Coulomb repulsion within the exter¬ 
nal harmonic trap, which orders the ions along the x-axis 
[351 [36] , and whose axial component reads 
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with uj the trap frequency along x. In the absence of the 
cavity the equilibrium positions {xf^^} form a chain, the 
interparticle distance dj = — x^^^ is inhomogeneous 

but almost uniform at the chain center [27], where it 
takes the minimal value d. Within the resonator the ions 
equilibrium positions {xj} are the minima of the total 
potential V = Idon + Kav- We choose the trap frequency 
UJ to ensure an incommensurate ratio between A and d, 
such that the dynamics are intrinsically frustrated. 

The strength of the cavity-mediated interactions is 
controlled by the cooperativity C = NU^/n^ which scales 
the strength of the nonlinear shift in Eq. Q. Eor |(7| ^ 1 
the mean photon number is independent of the ions posi¬ 
tions. In this limit Ecav ~ Pb cos‘^(kxj) and the total 
potential V can be mapped to the EK model [T21 EH [22] . 
The sliding-to-pinned transition is then expected at a 
critical value of the potential amplitude Vq oc and 
occurs at smaller values of rj with increasing |(7|, as il¬ 
lustrated in Eig. [^c). The sign of Uq, and thus of (7, 
determines the features of the pinned phase: Eor C > 0 
a pinned configuration minimizes Hat, since the minima 
of Pcav are at the nodes, while if (7 < 0 it maximizes Hat- 
As I(71 increases, the cavity potential changes shape, as 
in Fig. [^b), and it strongly depends on the value of B^ 
through the shift in Eq. §: Eor a fixed detuning Ac, the 
resonance Aeff = 0 is fulfilled for certain values of Hat, 
and thus for specific phases. Eor |Ac| > the resonance 
can directly separate the regime where the minima are 
either spikes or fiat bottomed. 

In order to evaluate the ions phase we define an ap¬ 
propriate thermodynamic limit: Since (7 oc A, we scale 
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FIG. 2 : Phase diagram as a function of C and 77 for (a) Ac = 0, (b) Ac = —2n. In (c) the phase diagram is plotted as a 
function of Ac and 77 for cooperativity C — — 2 . Parameters 77 and Ac are in units of The solid red line indicates the 
symmetry breaking transition from the sliding (S) to the pinned (P) phase, the hatched white region is an area of bistability. 
The color code gives the value of Bn for C > 0, and of 1 — Bn when C < 0. Subplots (d, e) display BnN vs. N for C — 0.5, 
Ac = 0, and (d) 77 = 50/^ and (e) 77 = 500/^ respectively, (f) The individual contributions cos^(kxj) to Bat for iV = 81 particles 
and deep in the pinned phase (Bn — 1.5 x 10“^) for C = 0.5. The ions are the parameters are |Ao| = 27r x 12 GHz, 

tv — 27t X 0.2MHz, for A = 11 the trap frequency is cj = 27r x 1.12MHz, while g is varied in order to sweep over different values 
of C. The cavity wavelength is A = 369nm and the ratio 2d/\ — 7.3507. The centre of the harmonic trap corresponds with a 
maximum of ITav, ensuring a symmetry breaking transition. The grey area in (a-c) indicates where the mean phonon number 
77 < 1 . Outside of this region our semiclassical analysis is reliable. 


Uq ^ 1/N to ensure that the resonance is at the same 
value of Bn [38]. To fix the ratio d/A, the trap frequency 
is scaled as cj ^ ^y\og(N)/N [36]. The equilibrium po¬ 
sitions of the ions in the total potential are numerically 
determined as a function of 77 , C, and A^ and their sta¬ 
bility is checked by means of a linear stability analysis. 
We characterize the stationary state by determining the 
corresponding bunching parameter and by classifying it 
in terms of sliding or pinned phase. For this purpose, we 
compare the minimum phonon frequency [22j [23] and an 
analog of the depinning force in the classical FK model 
iniESlEl], which is evaluated following the procedure 
discussed in Ref. [SUES- For C > 0, for instance, the 
trap center {x = 0 ) is at a maximum, so for an odd 
number of ions, the sliding to pinned transition can be 
observed through the sudden displacement of the central 
ion from the origin. For larger C a finite force is required 
to restore the symmetry of the system. At this transition 
we verify that the phonon gap is zero, and in the pinned 
regime it increases monotonically. Analogous considera¬ 
tions apply for (7 < 0 , where the maxima of the potential 
are at the antinodes of the field (in this case the center 
of the trap is shifted by A/4 to obtain a symmetry break¬ 
ing). We check that the values of 77 and C found by the 
symmetry breaking transition are the same at which the 
phonon gap starts to increase [32] , 

The resulting phase diagram is shown in Figs. Ha) 
and (b) for 11 ions and as a function of 77 and C for Ac = 
0 and Ac = — 2 /^ respectively, using the parameters of 
Ref. [ 12 ]. It exhibits sliding (S), pinned (P), and bistable 
phases (hatched). The color code gives the corresponding 
value of We have checked that the diagram remains 
substantially unvaried as the number of ions is scaled up 
according to the thermodynamic limit, apart from the 
bistable phases at (7 > 0, as we discuss below. We first 


consider the transition line, delimiting the S phase. This 
moves to smaller values of 77 as \C\ increases. For \C\ < 1 , 
the back-action due to the resonator can be neglected 
and the line follows the expected behaviour for the FK 
limit, with 77 = 77 c ^ Here, it separates the S 

from the P phase and exhibits the typical features of the 
Aubry transition. At larger values it changes functional 
dependence. Moreover, sliding and pinned phases can 
coexist for \C\ > 1 about the transition line, a typical 
feature of a first order transition. 

The bistable areas at (7 > 0, however, are of different 
nature than the one for (7 <0. For (7 > 0 they are due 
to finite-size effects. As C increases, in fact, the effective 
cavity potential becomes fiat except for the nodes, where 
it exhibits tight minima, see Fig. [^b). Thus, this poten¬ 
tial supports a stable sliding phase which is symmetric 
about the center, where the ions in general experience a 
fiat potential. As N is varied, bistability is still observed 
but it qualitatively changes its features. The bistable 
region in Fig. [^b) for (7 < 0 is instead due to the reso¬ 
nance Aeff = 0 and remains unvaried as N is scaled up 
according to our prescription. This resonance occurs for 
specific values of Bn^ and thus for specific sets of {%}. 
Analogous resonances have been reported in experiments 
with cold atoms in resonators jfiOHdT] and in theoretical 
works on similar setups [nnsnin]. In our case they indi¬ 
cate that either superlubric or stick-slip behaviors can be 
encountered depending on how the intensity 77 is varied in 
time. Fig. |^c) shows the phase diagram for C = —2 and 
as a function of 77 and Ac- For this parameter choice, the 
resonance Aeff = 0 exists only for with values 

varying between the minimum and the maximum value 
of — 2 Hat and thus of the bunching parameter. 

The bunching parameter is a particularly relevant 
quantity, since its value can be extracted from the in- 
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tensity of the light at the cavity output. In the sliding 
phase and for ^ 1, the particles are positioned at ev¬ 
ery point (modulus A) of the cavity potential and thus 
limTv^oo^AT = 0.5. In a commensurate phase ^ 0 
for C > 0 {Bn ^ 1 for C < 0 ). At small deviations 
from this limiting value, is a crude estimation of the 
kink density. The plots in Figs, [^a) and (b) show that 
Bn signals the transition from sliding to pinned when 
\C\ ^ 0.5, however this is not the case for larger C. 
So in general for Bn ^ 0.5 the phase is sliding, while 
for Bn < 0.05 at C > 0 (1 - Bat < 0.05 at C < 0) 
the phase is tightly pinned. The estimated kink number 
Nk ^ BnN grows linearly with as shown in Fig.j^d)- 
(e): The slope decreases as r] and C are increased, but 
never vanishes, thus showing that the pinned phase re¬ 
mains incommensurate. Nonetheless, deep in the pinned 
phase the ions at the chain edges, where the density is 
smaller, organize at commensurate distances with A, as 
depicted in Fig. [^f). This effect results from the choice 
of harmonic confinement and it shows that the edge ions 
enforce a boundary condition restricting the central ones 
from becoming truly commensurate. 

Further features of the cavity long range interaction 
manifest in the ions vibrations, as they give rise to fluc¬ 
tuations in the potential, which in turn affects the ions 
motion. We analyze them in the linear regime and we de¬ 
note 5a and 5xj as the quantum fluctuations of the cavity 
annihilation operator and of the ion positions about the 
mean values a = y/fi and %, respectively, where n is 
given in Eq. § and Xj are the minima of V. We de¬ 
compose the ions’ displacement Sxj in the normal modes 
Qn = (^n + ^n)/v^ Calculated at zero order in dd, with 
bn the bosonic operator annihilating a chain phonon at 
frequency Cavity and phonons are coupled by the 
linearized Heisenberg-Langevin equations in the presence 
of noise due to cavity decay and to an external damping 
reservoir coupled with the motion [nisi] 

5a — (iAgff ti)5a i<^ T 5]^) \/2^din i (5) 

bn — (i^n T r n^bn {5a T 5a) ) T \/2r n ^in,n : (b) 

where Cn denotes the cavity coupling with mode n 
and F^ is the mode’s damping rate. The Langevin 
operators (in = din,5in,n have zero mean value and 

([ 6 n (^05 Cin(^^ 0 ]) = ~ The solutions are station¬ 

ary when the eigenvalues possess no positive real parts. 
For F^ = 0 the stability is determined by the cavity pa¬ 
rameters and is warranted when Agff < 0 (the stability 
diagrams for the plots in Fig.[^a,b) are in the supplemen¬ 
tary material [32]). In this regime, retardation processes 
in photon scattering cool the chain to effective temper¬ 
atures T which depend on the detuning Ag and on the 
bunching parameter Bn- Thus, the ions stationary state 
determines the temperature at which the chain is cooled. 
In turn, for a given Ag disparate regions in Fig. [^a-c) 
are generally at different temperatures since Bn and C 
vary. From Eq. (§ one sees that for Ag < 0 and (7 > 0 
the ions are always cooled. Cooling for (7 < 0 is found 



FIG. 3: Spectrum at the cavity output S{h') (in arbitrary 
units) for (7 < 0: (a) in the sliding phase with Bn = 1 — 0.45 
and (b) in the pinned phase with Bn — ^ — 0.05 and for 
C = —0.5 (black line) and (7 = — 2 (red line), for 11 ^^^Yb^ 
ions, Ac = 0, Fn = 0.1/^, and T = lOO/xA. The resonances 
correspond to vibrational eigenmodes coupling with the cavity 
field and change in the pinned phase as C is increased. The 
elastic peak at ujp (corresponding to z/ = 0) is not reported. 


by suitably modifying the cavity detuning. In order to 
estimate T we cast the modes in the form of a covariance 
matrix: For the parameters of Ref. m and N = 11 ions 
we find that T ^ 125/iK can be reached. In this limit 
the standard deviation from the equilibrium positions, 
is smaller than A/2 in the pinned phase, thus 
the classical equilibrium positions dictate the phases of 
the system. This condition can be achieved for any point 
of the phase diagram, i.e., also for Agff > 0, by sym¬ 
pathetically cooling the chain [43], corresponding to an 
external reservoir with F^ > 0 in Eqs. (|2^-(23). 


The spectrum of the field emitted by the cavity, dout, 
contains information about the collective vibrational 
modes of the ions within the resonator im [45] . The out¬ 
put field is formally connected to the cavity field via the 
relation dout = din + v^^d [46] , and its power spectrum 
is given by S{v) cx {aont{v)'<hontiv)), where aout(i^) is the 
Fourier transform of (lout (t) [31j |32]. Figure displays 
S{v) for parameters such that the ions phases are sliding 
(a) and pinned (b). For each phase we took the same 
values of ^at but different values of (7. The peaks corre¬ 
spond to vibrational modes coupled to the cavity and it is 
apparent that in the pinned phase more peaks are visible. 
This is a result of the broken symmetry induced by the 
optical lattice potential. The effect of the cavity back- 
action is weak in the sliding phase as the only discernible 
change in the spectrum is its relative intensity. In the 
pinned phase, however, the intricacies of the back-action 
are particularly apparent. Here the spread of the cavity 
frequencies become more separated for (7 = —2 due to 
the softening of the cavity pinning (see Fig. [^b)), result¬ 
ing in the emergence of three distinct frequency bands. 
Contrary to this, when C = —0.5 the ions are tightly 
restricted to the potential minima resulting in a narrow 
frequency band. 

Our analysis is performed for parameters which are 
consistent with ongoing experiments, joining trapped 
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ions and CQED setups inHis] where the nonlinearity 
can be experimentally tuned by changing the number 
of atoms. This study is an example of competing long- 
range self-organization mechanisms which realizes a new 
paradigm of the Frenkel-Kontorova model. As well as dis¬ 
playing novel phases, such as bistability induced by the 
cavity mediated interactions, the inherent losses from the 
cavity can cool the ions in a controlled manner and al¬ 
low one to monitor the phases at the cavity output, thus 
setting the basis for feedback mechanisms controlling the 
thermodynamics of friction. 
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I. SUPPLEMENTARY INFORMATION 
A. The model 



FIG. 4: (a) The phonon gap plotted as a function of C with 
r] = IOOa^ (blue line) and r] = 150a^ (orange line). The vanish¬ 
ing of the phonon gap happens at the sliding to pinned tran¬ 
sition. (b) The phonon gap as a function of tj/k for C = 2.4 
showing bistable states where a sliding (black line) and pinned 
(red dashed line) phase co-exist. The parameters are the same 
as those in Fig.2 in the main article. 

and atomic variables are sufficiently small to justify the 
treatment. To this aim, we decompose the operators as 
a sum of mean values and fluctuations according to the 
prescription 

a = a 6a , 

Xj = Xj + 6xj , (9) 

Pj=Pi + ^pj > 

where (d) = d, {xj) = Xj, and {pj) = pj^ while the ex¬ 
pectation value of the fluctuations 6a,6xj,6pj vanishes. 


We treat the cavity mode in second quantization, with 
h = a the photon number operator, a the bosonic op¬ 
erator destroying a quantum at energy huJc and a) its 
adjoint. The resonator is driven by a classical field at 
frequency cjp, and we report the dynamics in the refer¬ 
ence frame rotating with The ions behave as classi¬ 
cally polarizable particles, which is valid when the atomic 
transition which couples with the cavity mode is far de¬ 
tuned from the fields, so that the detuning Aq = ujp — uJe\ 
is the largest frequency of the problem. Note that since 
|Ao| ^ |^c|: in the parameter regimes we consider, 
Ao ^ cjc — ujei. The coherent dynamics of N ions and 
cavity field mode is governed by the Hamilton operator 
{h=i) 


Mean values. The mean values satisfy the equations of 
motion 

9 

—a = (lAeff - Kja + T], 

d__pj 
m ’ 
d 

—pj = -^jFions - Uondj cos^(kxj ), 

with n = |dp and dj = djdxj the gradient with respect 
to the spatial coordinates of the j-th particle (evaluated 
at the equilibrium positions Xi,..., xw), while 


( 10 ) 

( 11 ) 

( 12 ) 


Aeff = Ac — UqNBn . (13) 


_^ p. _^ 

H = — cos^ kxj—i(pa—p'^a^). 

^ 2m 

j 3 

(7) 

The effect of cavity losses and damping on the motion is 
introduced by means of Heisenberg-Langevin (HL) equa¬ 
tion. For the moment we assume that the only incoherent 
effects are losses of the cavity field at rate so that the 
HL equation for the cavity field reads 

d = ^[d, i^] —/^d + \/^din(t). (8) 


In order to determine the classical equilibrium values we 
require that the quantities d, Xj and pj correspond to 
stationary solutions of the dynamical equations, namely 
dtCL = 0, dtXj = 0, and dtPj = 0. From ( p!Q| ) we obtain 


a = 


V 

K - iAeff 


(14) 


and with no loss of generality, we choose the phase of p 
such that a is real. Setting Eq. (11) to zero gives pj = 0. 
Substituting the value of n = |d|Mnto Eq. (12), one finds 
that it can be cast in the form 


We perform the study in the semiclassical limit, assum¬ 
ing that the fluctuations about the mean values of field 


— ^jAions djVcE 


(15) 
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with 

Uosm{kxj)cos{kxj) 

Oj i/cav ^2 1 + (A, - Uo Ee cos^kxeW/K^ ^ 

which gives 


where the derivatives in the expressions above are evalu¬ 
ated at the equilibrium positions. 

For convenience we introduce the normal modes of the 
crystal, that characterize the dynamics of the ions when 
the coupling with the quantum fluctuations of the cavity 
field can be neglected: 


Ir/P 

V^av = -arctan 



(17) 


n 




h . 

Qn •) 

mujn 


( 21 ) 


The equilibrium positions of the ions are then found by 
minimizing the total potential V = Idon + Kav such that 
the forces due to the Coulomb repulsion, the confining 
potential and the cavity field are balanced. 

Phonon gap. We determine the phonon gap by evalu¬ 
ating the dispersion relation of the ions vibrations in the 
total potential V (thus discarding the fluctuations of the 
cavity field but taking into account its mean effect). The 
phonon gap is plotted in Fig. |^a) for = 11 ions as a 
function of C. As noted in Ref. [ 22 ] the phonon gap does 
not vanish in the sliding phase as in the FK model due to 
the confining harmonic potential, however as the depth 
of the cavity potential is increased the phonon gap can 
fluctuate in the sliding phase. The transition from sliding 
to pinned can be visualized by the phonon gap vanishing 
at a critical cavity depth, after which the gap increases 
monotonically [29]. It should be noted that for an even 
number of ions a similar behaviour is observed however 
the phonon gap may not vanish around the transition 
point. This is a finite-size effect. 

The phonon gap across the bistable region is plotted 
in Fig. I^b) dependent on 77 in units of k. The black 
line shows an extended sliding phase due to the ions in 
the chain becoming stable as a result of the cavity back- 
action flattening the shape of the cavity potential. At the 
critical point the phonon gap increases abruptly signify¬ 
ing a drastic jump from the sliding to the pinned phase. 
The red dashed line shows the other stable configuration 
for the same paramaters, however in this case the sliding 
phase is shorter as the ions move to the pinned phase at 
a lower critical value of the cavity depth. 

Fluctuations. The coupled dynamics of the quantum 
fluctuations of field and motion are governed by the HL 
equations [46], which are found substituting the decom¬ 
position Eq. into Eq. and into the Heisenberg 
equations of motion for the center-of-mass variables, and 
using that the mean values are the stationary solutions. 
The equations read 

Sa = (iAeff — K)Sa — iUoa cos‘^{kxi) + \/^din , 

i 

(18) 

Sxj = , (19) 

m 

6pj = - ^{6xidi) {djVions + Uondj cos‘^{kxj)) 
i 

— Uq (d*M + a6a^) dj cos^{kxj ), (20) 


with Mjn the element of the orthogonal matrix relating 
the local coordinates 6xj with the normal-mode coordi¬ 
nates, that diagonalize Eqs. ([T^-( 2 Q) when the cavity 
fluctuations 6a are set to zero; the variables are the 
dimensionless position coordinates of the normal modes. 

We denote by bn and the bosonic operators annihi¬ 
lating and creating, respectively, a phonon of the normal 
mode at frequency They are defined through the 
equations = ( 6 „ + bl)/\/2 and = i(bl - bn)/V^, 
and the dynamical equations for the fluctuations then 
take the form: 


Sa = (iAeff - i^)Sa - la'^CniK + bl^) + V^din , (22) 

n 

bn — (i^n T ^iaCn(^Sa T ^d^) “I" 2r n bin,n 5 

(23) 

which also includes the coupling of mode n to a reser¬ 
voir at rate F^. The corresponding Langevin force is de¬ 
scribed by the input noise operator bi^^n^ with ( 6 in,n) = 0 
and 


nit') kn,n'{t")) = “ t") ■ (24) 


with Nn = N{ujn) the mean excitation number of an 
oscillator of frequency ujn at the temperature of the con¬ 
sidered environment. The coefficients Cn in Eq. (22)-(23) 
read 


Cn = ]J 2mto E cos^ikxj), (25) 

where the derivatives are evaluated at the equilibrium 
positions Xj. 

Stability diagrams and cavity cooling The stability of 
the system is dependent on the equations of motion cou¬ 
pling the cavity and motional fluctuations which can be 
written in the compact form 

^ + tut) (26) 

where the elements of the matrix A contain the argu¬ 
ments of Eq. ( [^ and Eq. ( [^ [31]. Eor the system to 
be stable the real parts of the eigenvalues of A must be 
negative, otherwise at least one mode of the chain will be 
heated. It can be shown that the stability condition is 
Aeff < 0. The stable regions of our parameter space are 
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200 400 600 800 200 400 600 800 

tjIk, rf/hi 

FIG. 5: Unstable (yellow) and stable (blue) regions of the 
system for (a) Ac = 0 and (b) Ac = — 2/^. For the system to 
be stable the eigenvalues of the matrix A (see Eq. 42 in M) 
must have negative real parts. If this is not the case then 
there will be at least one mode which is heated by the cavity. 



200 400 600 800 500 1500 2500 

r]/K r]/K 


FIG. 6 : Average chain temperature in units of degrees Kelvin 
versus C and 77 in units of kappa. The chain is cooled in 
the colored area by the cavity, which is reached by suitably 
modifying the detuning, shown for (a) Ac = —2k and (b) 
Ac = — IOa^. The solid red line is the symmetry breaking 
transition, and the white area is the region of the parameter 
space which experiences heating from the cavity. 

plotted in Fig. for zero and finite A^, with F^ = OVn. 
In the extended stable region produced by the detuning 
the chain is cooled by the cavity to sufficiently low tem¬ 
peratures (minimum of around 125/iK for Ac = — IO/ 5 :) 
and this is shown in Figj^ The cavity cooled region may 
be increased by selecting a larger Ac, and the transition 
may be observed as a function of 77 by choosing specific 
C thereby only experiencing a small increase in temper¬ 
ature. 

In order to stabilise the rest of the unstable region, 
all the modes can be coupled to an environment that 
can damp the excitations (this coupling is denoted by 
a finite F^). To stabilize the entirety of the unstable 
regions shown here we must choose the coupling to be 
F ^ = 0.1k with the temperature of the external bath 
being Text = lOOjaK. 


xlO'"* 



FIG. 7: Restoring force in units of muj‘^L versus C, where 
L — and is of the order of the interpar¬ 

ticle distance at the chain center d. Three different pump 
strengths are shown (from left to right): 77 = IOO/ 5 : (blue), 
77 = 105/5: (red) and 77 = 110k (yellow) (the critical value of G 
becomes smaller with increasing 77 ). The symmetry breaking 
transition is also shown as a black dashed line for each case. 
The parameters are the same as those in Fig.2 in the main 
article. 

to the de-pinning force of the Frenkel-Kontorova model, 
whereby in the sliding phase with the symmetry of the 
system intact the restoring force is zero. However in the 
pinned phase the restoring force becomes finite thereby 
realizing the growth of static friction in the system. A 
finite force is applied to all the ions equally which may 
be physically implemented by tilting the lattice potential 
[24]. The sum of the total forces acting on all the ions 
must vanish to ensure a stable solution and are given by 

F = Tion + Fcav + Ft = 0 (27) 

where Fio„ = and Fcav = and Ft is the 

finite force applied to the chain. The critical value of the 
restoring force is then taken as the finite force required to 
move the central ion back to the maximum of the cavity 
field, i.e. such that cos^(/cxo) = 1 at which point Tres = 
Ft- The restoring force is shown in Fig.j^as a function of 
C for different pump strengths p/k. Before the transition 
the restoring force is zero, however after the transition 
the restoring force becomes finite highlighting the broken 
symmetry of the system. 


B. Evaluation of the restoring force 

C. Spectrum at the cavity output 

To discern the sliding and pinned phases we calculate 

the restoring force required to return the central ion back We report here the analytical results, which were ob- 
to the maxima of the cavity field. This is analogous tained in Ref. m- The spectrum is calculated after 























evaluating the Fourier transform of the linearized HL for 
the fluctuations, using a = d + ^a. The quantum compo¬ 
nent of the spectrum reads 


Of,.) _ 


(28) 


where the first term is due to coupling of the quan¬ 
tum vacuum with the crystal vibrations with d(z/) = 

2 

2 I • ^2 and the second term is due to thermal 
noise coupling to the modes. The prefactor is given by 


where Sd(iy) is the Fourier transform of Sd, Sd(iy)^ is the 
Hermitian conjugate of Sd(iy), and we omit the Rayleigh 
peak at z/ = 0, i.e., uj = ujp^ which corresponds to the 
classical part. After some algebra, one finds the expres¬ 
sions 


s{p) = S'o(i 2 )| 


4K\e{p}\^a^ 

+ iiy- Aeff)2 

2r2 ^ Vi 




+ 


Nn + 1 


Tl + {u:n-vY Tl + iujn + i'r 


So{v) = 


+ (jy + Aeff)2 


1 + 


46»(i/)Aeffa' 


7^2 


{k - ivY + A^ 


eff 


(30) 

(29) These are the formula we employed to evaluate Fig. 3. 
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